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It is argued that various kind of shell structure which oc- 
curs at specific values of the magnetic field should be observ- 
able in small quantum dots in the addition energies and the 
magnetic susceptibility. 

PACS number(s):73.20Dx, 73.23Ps 



Quantum dots are ideal mesoscopic objects for the 
study of quantum mechanical properties and comparison 
with classical behavior |l| . The smaller the size of quan- 
tum dots, the larger the prevalence of quantum effects 
upon the static and dynamic properties of these systems. 
Quantum dots may be considered as true artificial atoms 
the properties of which can be controlled by men. 

An important feature of finite Fermi systems is the ex- 
istence of shells which give rise to magic numbers. The 
periodic table of chemical elements is a striking exam- 
ple for the dramatic difference between closed shell and 
partially filled shell systems; for nuclei and metallic clus- 
ters see for example |^-^. In recent experiments |^|^] 
shell structure phenomena have been observed clearly for 
quantum dots. In particular, the energy needed to place 
the extra electron (addition energy) into a vertical quan- 
tum dot has characteristic maxima which correspond to 
the sequence of magic numbers of a two-dimensional har- 
monic oscillator In this Letter we propose manifes- 
tations of shell effects that should be observable in the 
addition energy and the magnetic susceptibility for small 
quantum dots under the influence of a magnetic field. 

We choose as the confining mean field for the electrons 
in quantum dots the harmonic oscillator potential. Any 
smooth and finite potential admitting bound states can, 
for the lowest few levels, generically be approximated by 
the harmonic oscillator potential For small electron 
numbers the harmonic oscillator has been used success- 
fully as a phenomenological effective confining potential 
in quantum dots As was shown in the simple case of 
a two-electron system, the external parabolic potential 
and the Coulomb interaction leads to a new oscillator 
frequency of the effective mean field |jT^. The effect of 
an external homogeneous magnetic field can be calcu- 
lated exactly for a three dimensional (3D) harmonic os- 
cillator potential irrespective of the direction of the field 
pl| . Our discussion here is based upon the 2D version 
of the Hamiltonian [|ll[ . The magnetic field acts perpen- 
dicular to the plane of motion of the electrons and the 
spin degree of freedom is incorporated, i.e. H = X] i=i 



where 



h = ^{P- -Af + '^{u;lx^ + u^y) + ^V.S.. (I) 
2m* > a 
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Here A = [r x B]/2 and az is the Pauli matrix. 
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FIG. I. Single-particle spectra as a function of the mag- 
netic field strength. Spectra are displayed for (a) a plain 
isotropic two dimensional oscillator, (b) a deformed and (c) 
an isotropic oscillator including an L^-term. For better illus- 
tration the value 2/i* is used for the spin magnetic splitting 
in all Figures. 

We do not take into account the effect of finite tem- 
perature; this is appropriate for experiments which are 
performed at temperatures kT <^ A with A being the 
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mean level spacing. In the following we use meV for the 
energy and Tesla for the magnetic field strength. 

In Fig. la we display the single-particle spectrum for 
the isotropic case {ux = ^y) as a function of the mag- 
netic field strength B . The effective mass which deter- 
mines the orbital magnetic moment for the electrons is 
chosen as m* — 0.067TOe. It leads to /i'^ — 15/is while 
the effective spin magnetic moment is fj,* — 0.5fj,B- The 
magnetic orbital effect is much enhanced in comparison 
with the magnetic spin effect, yet the tiny spin splitting 
does produce signatures as we see below. 

Shell structure occurs whenever the ratio of the two 
eigenmodes V,± of the Hamiltonian (|^) (see Ref. is 
a rational number with a small numerator and denom- 
inator. The shell structure is particularly pronounced 
if the ratio is equal to one (for B = and ujx — ^y) 
or two (for B « 1.23 ) or three (for B « 2.01 ) and 
lesser pronounced if the ratio is 3/2 (for B = 0.72) or 
5/2 (for B = 1.65). The values given here for B depend 



on m* and 



As a consequence, a material with an 



even smaller effective mass m* would show these effects 
for a correspondingly smaller magnetic field. The magic 
numbers (including spin) turn out, for _B = 0, to be the 
usual sequence of the two dimensional isotropic oscilla- 
tor, that is 2,6,12,20,... §]. For B « 1.23 we find a 
new shell structure as if the confining potential would be 
a deformed harmonic oscillator without magnetic field. 
The magic numbers are 2,4,8,12,18,24,... which are 
just the numbers obtained from the two dimensional os- 
cillator with Lu^ = 2lu^ (cli> and a;< denote the larger and 
smaller value of the two frequencies). Similarly, we get 
for B « 2.01 the magic numbers 2, 4, 6, 10, 14, 18, 24, . . . 
which corresponds to cjy = 3(jJ<. 

If we start from the outset with a deformed mean field, 
i.e. if we choose, say, uix = — P)^y with /3 > 
two major effects are found: (i) the degeneracies (shell 
structure) are lifted at i? = depending on the actual 
value of /3, and (ii) the values for B at which the new shell 
structures occur are shifted to lower values. In Fig. lb we 
display an example referring to (3 = 0.2. The significance 
of this finding lies in the restoration of shell structures by 
the magnetic field in an isolated quantum dot that does 
not give rise to magic numbers at zero field strength due 
to deformation. We mention that the choice /3 = 0.5 
would shift the pattern found at _B w 1.23 in Fig. la to 
the value B = 0. 

It is the shell structure caused by the effective mean 
field which produces the maxima that are observed ex- 
perimentally in the addition energy /i(iV -I- 1) — fi{N) — 
Ejs[+i — En + /C Here En is the single-particle en- 
ergy of the effective mean field in quantum dots, e^/C is 
the electrostatic energy and /i(iV) is the chemical poten- 
tial. The electrostatic energy is much larger than the 
difference En+i — En, however, it is the fluctuations 
(shell effects) of the difference that matters, at least for 
small quantum dots Quantized electronic energies 



have been observed even for nanometer-scale metal par- 
ticles JT5| . A similar effect is known in nuclear physics 
and for metallic clusters. There, shell effects due to 
single-particle motion create minima in the total poten- 
tial energy surface which is dominated by the bulk energy, 
which is the classical liquid drop energy ||-|^ . The anal- 
ogy goes further in that, in an isolated small quantum 
dot, the external magnetic field acts like the rotation on 
a nucleus thus creating new shell structure; in this way 
superdeformation (axis ratio 2:1) has been established 
for rotating nuclei owing to the shell gaps in the single- 
particle spectrum ||]. 

In a previous study ||ll[ we have obtained various 
shapes of the quantum dot by energy minimization. In 
this context it is worth noting that at the particular val- 
ues of the magnetic field, where a pronounced shell struc- 
ture occurs, the energy minimum would be obtained for 
circular dots, if the particle number is chosen to be equal 
to the magic numbers. Deviations from those magic num- 
bers usually give rise to deformed shapes at the energy 
minimum. To what extent these 'spontaneous' defor- 
mations actually occur (which is the Jahn-Teller effect 
Jl6t), is subject to more detailed experimental informa- 
tion. The far-infared spectroscopy in a small isolated 
quantum dot could be a useful tool to provide pertinent 
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FIG. 2. Magnetic susceptibility x ~ —d^Etot/dB^ in arbi- 
trary units as a function of the magnetic field strength for the 
isotropic oscillator without L^-term. 

The question arises as to what extent our findings de- 
pend upon the particular choice of the mean field. Here 
we confirm the qualitative argument presented above 
that for sufficiently low electron numbers virtually any 
binding potential will produce the patterns found for the 
harmonic oscillator. We add to the Hamiltonian (|l|) the 
term —XhuL^ where L is the dimensionless z-component 
of the angular momentum operator. In nuclear physics, 
the three dimensional analogue of the combined prob- 
lem is known as the Nilsson model and has been quite 
successful in explaining the spectra of deformed nuclei 
Q . For A > the additional term lowers the energy lev- 
els of higher angular momenta, in other words, its effect 
mimics a bulging out of the lower part of the harmonic 
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oscillator potential. As a consequence, it has the effect of 
interpolating between the oscillator and the square well 
single-particle spectrum ||]. For uj^ ^ uiy and A 7^ the 
combined Hamiltonian H' = H—XhcoL^ is non-integrable 
I p^ and the level crossings encountered in Figs.l become 
avoided level crossings. The essential effect upon the 
lower end of the spectrum can be seen in the isotropic os- 
cillator where the magnetic quantum number m is a good 
quantum number. In this case H' = H^^°^'' — XfiLom? . 

In Fig.lc we display a spectrum of such H' . The shell 
structure, which prevails for A = throughout the spec- 
trum at _B « 1.23 or B K, 2.01, is now disturbed to an 
increasing extent with increasing shell number. However, 
for the parameters chosen the structure is still clearly dis- 
cernible for about seven shells, that is for particle num- 
bers up to about twenty five. The lifting of the degenera- 
cies at B = is also clearly seen where the levels are split 
according to the absolute values of |m|; it is this splitting 
which gives us guidance in choosing an appropriate value 
for A: for i? = we aim at levels that lie between the cor- 
responding degenerate levels pertaining to the harmonic 
oscillator and the two dimensional square well where the 
splitting of these levels is very strong. 

(a) 

% 

N=14 

[ 

4 
2 



N=15 

[ 

4 
2 



quantum dot of fixed electron number, the ground state 
will undergo a rearrangement at the values of where 
level crossings occur. This should be observable in the 
magnetic susceptibility as it is proportional to the second 
derivative of the total energy with respect to the field 
strength. In Figs. 2 and 3 we have plotted —d^Etot/dB'^ 
where i?tot is the sum of the single-particle energies filled 
from the bottom up to N which is the electron number. 
We discern clearly distinct patterns depending on the 
electron number, in fact, the susceptibility appears to be 
a fingerprint pertaining to the electron number. Fig. 2 
and Fig. 3a illustrate results for A = and Fig. 3b for 
A = 0.01. The deformed oscillator does not produce new 
features except for the fact that all lines in Figs. 2 and 3 
would be shifted towards the left in accordance with the 
discussion above. 

All features of Figs. 2 and 3 can be understood from the 
single-particle spectra displayed in Fig. la and Fig.lc. If 
there is no level crossing, the second derivative of i?tot is 
a smooth function. The crossing of two occupied levels 
does not change the smoothness. In contrast, if the last 
occupied level takes part in a level crossing in such a 
way that it is ascending before and descending after the 
crossing, the second derivative of i?tot must show a spike. 
In other words, spikes occur if the last occupied level is 
bent concavely when viewed from the abscissa; a convex 
curvature has no effect. As a consequence, we understand 
the even-odd effect when comparing N = 8 with = 9 
in Fig.2 and iV = 14 with A^ = 15 in Fig. 3. The spin 
splitting caused by the magnetic field at B w 2.01 for 
A^ = 8 is absent for Af = 9. 
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FIG. 3. Similar to Fig.2 for different particle numbers and 
(a) for A = and (b) for A = 0.01. 

When the magnetic field is changed continuously for a 




FIG. 4. Blow ups of the relevant level crossings explaining 
(a) the features in Fig.2 and (b) those of Fig. 3. The left 
and right hand drawing of (a) refers to A'^ = 8 and = 9, 
respectively. The left hand and right hand drawing of (b) 
refers to A = and A = 0.01, respectively; here the thick and 
dotted lines refer to A'^ = 14 and A*' — 15. 

This becomes evident when looking at a blow up of 
this particular level crossing which is illustrated in Fig. 4a, 
where the last occupied level is indicated as a thick line 
and the points where a spike occurs are indicated by a 
dot. In this way, we also understand that a splitting 
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occurs for A'^ = 6 at _B w 1.23 (not presented here) and 
no spike at all for TV = 8 at that same field strength. 
Note that the splitting is proportional to the effective 
spin magnetic moment Another illustration of the 
even-odd effect, which persists through all numbers, is 
given in Fig. 3. Now we obtain a triple splitting for N = 
14 at B Ri 1.23 which becomes a double splitting for N — 
15, while the double splitting for N = lA at B 2.61 
becomes a single line for N — 15. The double splitting 
is associated with the crossing of only two double lines, 
while the triple splitting originates from the crossing of 
three double lines. The line distance at i? « 2.61 is 
increased as it occurs at a larger field strength. 

When we consider the crossing of only two double lines 
as is the case for = 8 and = 9, or = 14 and 
N = 15 at B Ki 2.61, the effect of the L^-term is unim- 
portant. However, the L^-term does become significant 
when we turn to level crossings of three or more dou- 
ble lines, in fact, the term removes the artificial pattern 
of level crossings of three and more lines. Wc focus our 
attention at the multiple crossings around B = 1.2 by 
comparing Fig. 3a with Fig. 3b. In Fig. 4b we illustrate a 
blow up of the situation which produces the different re- 
sults for A and A = 0.01. For A 0.01 the levels 
with a non-zero slope are lowered as they refer to m 7^ 
while the horizontal levels remain unaffected, hence the 
crossings change as seen by comparing the left and right 
hand drawing of Fig. 4b. The last occupied level is indi- 
cated by a thick line and a dotted line for the A^ = 14 
and A^ = 15 system, respectively; the corresponding po- 
sitions, where spikes occur, are marked by solid dots and 
open circles. The even-odd effect is encountered again 
but appears converted for A — 0.01. Note, however, that 
the triple line for A^ = 15 at i? = 1.2 is not a spin split- 
ting but an effect due to the L^-term; in other words, 
for /i* ^ the double line distance of A^ = 14 in Fig. 3b 
would vanish while the triple line would not as it depends 
on A. Spin flips can also be understood by the same to- 
ken, here the general rule applies: the crossing of the 
top (bottom) with the bottom (top) line of a double line 
causes the spin to flip. Hence, both lines of the double 
splitting in Fig. 2 are associated with a spin flip (A^ — 8), 
but neither of the single lines (A^ — 9). Also, in Fig. 3b 
the triplet (A^ — 15) has no spin flip but if A would be 
chosen sufficiently small, the central line of the triplet 
would be split further and a spin flip would be associ- 
ated with the doublet (which is then part of a quartet). 
Strictly speaking, the spikes are (5-f unctions with a factor 
which is determined by the angle at which the two rele- 
vant lines cross. Our flgures are numerical results which 
do not exactly reflect this feature. If the level crossings 
are replaced by avoided crossings (Landau-Zener cross- 
ings), the lines would be broadened. This would be the 
case in the present model for A > and (3 > 0. Finite 
temperature will also result in line broadening. 

To summarize: the consequences of shell structure ef- 



fects for the addition energy of a small isolated quantum 
dot have been analyzed. At certain values of the mag- 
netic field strength shell structures appear in a quan- 
tum dot, also in cases where deformation does not give 
rise to magic numbers at zero field strength. Measure- 
ments of the magnetic susceptibility are expected to re- 
flect the properties of the single-particle spectrum and 
should display characteristic patterns depending on the 
particle number. The latter property could be of interest 
in applications because it enables control of the electron 
number in small isolated quantum dots and the magnetic 
held strength. The splittings due to the spin magnetic 
moment provide for a quantitative assessment of the ef- 
fective spin magnetic moment fi* . 
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